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ABSTRACT 



Since gravitational lensing effects directly probe inhomogeneities of dark mat- 
ter, lensing-galaxy cross- correlations can provide us important information on the 
relation between dark matter and galaxy distributions, i.e., the bias. In this pa- 
per, we propose a method to measure the stochasticity/nonlinearity of the galaxy 
bias through correlation studies of the cosmic shear and galaxy number fluctu- 
ations. Specifically, we employ the aperture mass statistics Map to describe the 
cosmic shear. We divide the foreground galaxy redshift Zf < Zg into several 
bins, where Zg is the redshift of the source galaxies, and calculate the quantity 
< MapNg{zf) >^ / < Ng{zf) > for each redshift bin. Then the ratio of the 
summation of < MapNg{zf) >^ / < Ng{zf) > over the bins to < Mjp > gives a 
measure of the nonlinear/stochastic bias. Here Ng{zf) is the projected surface 
number density fluctuation of foreground galaxies at redshift Zf, and Map is the 
aperture mass from the cosmic-shear analysis. We estimate that for a moder- 
ately deep weak-lensing survey with Zg = I, source galaxy surface number density 
rii) = 30 gal/ arcmin^ and a survey area of 25 deg^, the effective r-parameter that 
represents the deviation from the linear and deterministic bias is detectable in 
the angular range of I'-IO' if |r — 1| > 10%. For shallow, wide surveys such as 
the Sloan Digital Sky Survey with Zg = 0.5, = 5 gal/ arcmin^, and a survey 
area of 10^ deg^ a 10% detection of r is possible over the angular range 1' — 100'. 

Subject headings: cosmology: theory — dark matter — galaxy: cluster: general — 
gravitational lensing — large-scale structure of universe 
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1. Introduction 

Galaxy surveys have provided us abundant knowledge of the large-scale structures of 
the universe. To extract cosmological information from galaxy surveys, however, it is crucial 
to understand the bias, i.e., the relation between the galaxy distribution and the underlying 
dark matter distribution (e.g., Strauss 1999 and references therein). The simplest bias model 
assumes a constant (in scale) bias factor between the two-point correlation of galaxies and 
that of dark matter: i.e., 

Ur)-b'arl (1) 

where ^gg is the two-point correlation function of galaxies and ^ is the two-point auto- 
correlation of dark matter. The 6-factor is referred to as the bias factor, which can be a 
function of the galaxy type, brightness, etc.. The above relation is further extended to 

Sg = bS, (2) 

where Sg — Sn/n, the galaxy number density fluctuation, and 5 — Sp/ p, the dark matter mass 
density fluctuation. Note that relation (2) is more restrictive than relation (1) and is referred 
to as the linear and deterministic bias relation. Under this assumption, it is straightforward 
to get the dark matter clustering properties from those of galaxies since only the amplitude 
of clustering needs to be adjusted. In reality, the bias could be much more complex than this. 
The bias factor b can be scale-dependent, which would introduce complexities in converting 
the galaxy distribution to the underlying dark matter distribution. Further more, since 
galaxy formation is a complicated process, it is expected that there should be a certain level 
of nonlinearity and stochasticity in the relative distribution between the dark matter and 
galaxies (e.g., Dekel & Lahav 1999). Indeed the deviations from the hnear and deterministic 
bias have been seen in numerical simulations (e.g., Dekel & Lahav 1999; Blanton et al. 1999). 
It is shown that older galaxies and the underlying dark matter are highly correlated, but 
younger galaxies are poorly correlated with dark matter (e.g., Blanton et al. 1999, 2000; 
Somerville et al. 2001) 

As far as the second moments are concerned, the nonlinear/stochastic bias can be de- 
scribed by the quantities b and r which are defined as 

(7^ 



and 



^ f,, (3) 



iMi, (4) 
aga 



where =< > and =< S"^ >. In the case of linear and deterministic bias, r — 1. 
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It has been proposed to measure the nonhnearity/stochasticity of the bias using the 
information of redshift distortions of the galaxy distribution, which arc presumably caused 
by the underlying dark matter distribution (Dekel & Lahav 1999; Pen 1998). Tegmark 
and Bromley (1999) analyzed galaxy correlations in the Las Campanas Redshift Survey ( 
Shectman et al. 1996; Bromley et al. 1998). They found that different subpopulations of 
galaxies are not correlated perfectly with each other, indicating that the galaxy densities 
and the underlying dark matter density cannot be perfectly correlated, or in other words, 
the bias is nonlinear/stochastic. They further showed that the deviation from the linear and 
deterministic bias is significant for late-type galaxies, and the r factor ranges from ~ 0.98 
for very early-type galaxies to ~ 0.57 for late-type ones. However, Blanton (2000) pointed 
out that this large stochasticity (small r) could be the result of selection effects, and his 
estimation of r is r w 0.95. 

Gravitational lensing effects due to large-scale structures (cosmic shear) have been be- 
coming important probes for studying the dark matter distribution in the universe (e.g., 
Blandford et al. 1991; Kaiser 1992; Schneider et al. 1998; Melher 1999; Bartelmann & 
Schneider 2001). Although weak (a/< 7^ > is on the order of ~ 1% on angular scales < 10', 
where < 7^ > is the top-hat variance of the cosmic shear), the cosmic-shear signals have 
been detected by different groups (e.g., Mellier et al. 2002 and the references there in), and 
useful constraints on (jg and flo have been set up with these lensing analyses (e.g., Maoli 
et al. 2001). Recently, galaxy-galaxy lensing has also been detected (Fisher et al. 2000; 
McKay et al. 2001; Hoekstra et al. 2002), and this enables the statistical study of properties 
of galactic halos and galaxy-mass correlations (e.g, Guzik & Seljak 2001, 2002). 

While the most direct gain from cosmic-shear analyses is the projected dark matter 
distribution, or the projected power spectrum in Fourier domain (Bartelmann & Schneider 
1999; van Waerbeke et al. 1999), the cross-correlation between the cosmic shear and the 
galaxy distribution refiects the relative distribution of dark matter and galaxies, i.e., the 
bias. The combination of this cross-correlation and the galaxy autocorrelation has been 
used to measure the scale dependency of the 6-factor (van Waerbeke 1998; Hoekstra et al. 
2001). Hoekstra et al. (2002) also made a first try to measure the r-factor by using the cosmic 
shear-galaxy, galaxy-galaxy, and cosmic shear-cosmic shear correlations. In their sample, the 
foreground galaxies occupy a relatively narrow redshift range. On the other hand, the mass 
contribution to the cosmic shear comes from a relatively broad redshift range. Therefore 
they need a conversion factor to convert the direct measurement result to the r-factor (see 
eq.[17] in Hoekstra et al. 2002). The conversion factor is model-dependent, and this adds 
complications to the interpretation of the observational results. 

With fast developments in observations, it will be possible to divide galaxies according 
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to their redshifts, which can be obtained at least photometrically. This will enhance the 
power to extract cosmological information from lensing analyses. 

In this paper, we propose to use foreground galaxies in several redshift bins with Zf < Zg 
where Zs is the source redshift, to get a summation of < MapNg{zf) / < Ng{zf) > over the 
bins, which should give a good estimation of < M^^ > if the bias is linear and deterministic, 
i.e., r = 1. In other words, the ratio [S^^. < MapNg{zf) >^ / < Ng{zf) >]/ < M^^ > 
measures the projected r^-factor weighted by the power spectrum in different redshift bins. 
Here Map is the aperture mass used to describe cosmic shear signals, and Ng is the surface 
number density fluctuations of foreground galaxies. This measurement is largely independent 
of cosmology and of specific types of power spectra, and there is no conversion factor between 
r and the direct measurement result. 

The rest of the paper is organized as follows. Section 2 discusses the effects of different 
binning of redshift of foreground galaxies on the r-estimation, i.e., the intrinsic uncertainties. 
In section 3 we estimate the signal-to-noise level of the measurement of r for various models. 
Section 4 includes discussions. 



2. Redshift Binning 

In this paper, the aperture-mass statistics are used in the cosmic-shear analysis (e.g., 
Schneider et al. 1998). The aperture mass MapiOc) is defined as the projected surface mass 
density filtered with a compensated filter function (Schneider et al. 1998) 

Map{9c) = J d''9U{9,ec)K{9), (5) 

where k, is the dimensionless surface mass density, U is a compensated filter function with 
U — for 9 > 9c- The quantity M^p is related to the observable shear 7 through 

Map{9c) = J d^9Q{9,9cht{9), (6) 

where jt is the tangential component of the shear 7, and Q{9, 9c) = (2/6*^) d9'9'U{9' , 9c) — 
U{9,9c). 

The surface mass density k is associated with the density perturbation along the line of 
sight by 

3/i/o\2 ^ 5{fK{w)9,w) 

W = o — ] ^0 dwg{w)fK{w) , (7) 

2 \ c / Jo a{w) 
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where Hq is the Hubble constant, is the current matter density parameter of the universe, 
c is the speed of hght, w is the comoving radial distance, wh is the comoving radial distance 
to the horizon, fx is the comoving angular diameter distance, a is the scale factor, S is the 
density perturbation, and g is defined as 

g(w)= dwp^{w) — , (8) 

Jw JK{W') 

where Pw{w) is the distance distribution of source galaxies, and Pw{w)dw = Pz{z)dz, where 
Pz{z) is their corresponding redshift-distribution. We use the filter function of Schneider 
et al. (1998), which has the form C/(^, ^c) = u{e 10^)101 with u{x) = for x > 1, and 

and 

,M.(l±i)£±i),=(i_,^)', 

TT 

with / an integer. 

Then in the Fourier domain, we have 

where P is the three-dimensional power spectrum of density fluctuations at the time corre- 
sponding to w, and /; is defined as 

Ii{ri) = dx X u{x)Jo{xri) = -ri'^^+^Ui+3{ri). 

Jo TT 

In the following we will use / = 1, and then 

7(77) = '^V-'Mv)- (10) 



TT 



In comparison with the top-hat window function, the compensated filter function (10) is 
very localized with a peak at 77 ~ 4, and therefore < M^p{6c) > provides a good measurement 
on the power spectrum of projected density fiuctuations at scale s ~ (Bartelmann & 
Schneider 1999). When the cross-correlation between Map and the foreground-galaxy distri- 
bution is considered, the narrowness of the redshift distribution of foreground galaxies makes 
the cross-correlation signals come from a well-defined cosmic time. Thus, in conjunction with 
P, the contributions of the mass distribution to the cross-correlation peak at a well-defined 
scale k ^ {A/9c)/ fxizf), where Zf is the peak redshift of foreground galaxies and fxizf) is 
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the angular diameter distance to Zf. This property has been used in studying the scale de- 
pendence oir/h through the measurement of the lensing-foreground-galaxy cross-correlation 
and the auto-correlation of foreground-galaxy distributions (Schneider 1998, van Waerbcke 
1998). In the analysis presented in this paper, we also use this property as explained in the 
following. 

The projected distribution of foreground galaxies is described by 

(11) 

where A^(^^) is the surface number density of galaxies in the direction 9^ and N is the mean 
surface number density of galaxies. Then the cross-correlation of M^p and Ng is 

< M,p{e,)Ng{e,) >= 37r(^)'fio b r [ dw^^^^^f^-^ [ ds s pf^^^,w)r\se,), (12) 

V c / J a{w)fK{w) J ^Jk{w) ) 



and the auto correlation of Ng is 



< N]{Q,) >= 27r62 / dw^^-Jri—^ I ds s p(-^ w)p(se,), (13) 

J Jk\w) J \fK[w) J 

where P/ is the foreground-galaxy distribution, h and r are the bias factors discussed in the 
section 1, and r = 1 if the bias is linear and deterministic. 



As stated before, if p/ is highly peaked at Zf, then 

< MaM)Ng{ec) > oc 37r(^)'oo b r—^^^l^P{k,Zf), (14) 

\ c / a[Zf)fK[Zf) 

and 



< Ng\^c) > oc 27r6 -2"^- '-P{k,Zf), (15) 

JK\^f) 

with k {4:/9c)/ fxizf). Notice that in equations (14) and (15), we use Zf instead of w 
to label the time, and the formulation for g and fx should be changed accordingly in the 
numerical calculations. Then 



< N^iOc) > "^'^ I 2 \ c J "a''{zf) ^uzyzf 

Prom equation (9), it is seen that 

<M4(«j>^[l(f)'"5/<^.f;i(S/(*-) 
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with k {4:/9c)/ fxiz). Therefore if the foreground galaxies are divided into several redshift 
bins with Zf < Zg, and [< Map{0c)Ng{9c) >]^/ < Ng{9c) > is calculated for each bin, the 
summation of [< Map{9c)Ng{9c) >]^/ < N'^{9c) > xp/(-2/) over the bins would give rise to 
an estimation of < M'^^{9c) > if r = 1. With r 7^ 1, then we have 



E 



^ [<M,p{9,)Ng{9,,Zf) >J- \^r^^ M^(9\^ 



< m(9,, Zf) > 



or 



^ E,,{[< M,MNg{9,,Zf) >]V < N^g{9,,zf) > xpf{zf)Azf} 

Here we have explicitly written down the redshift dependence of Ng. The approximate sign " 
~ " means that when equation (17) is used in estimating r, there are intrinsic uncertainties 
depending on the redshift binning. The binning effects are discussed below. In equation 
(17), it is assumed that r is independent of scales. In general, r is scale-dependent, and then 
equation (17) provides a measurement of the power-spectrum- weighted r- factor. In addition, 
theoretical studies and numerical simulations showed that the r factor evolves with time, 
i.e., it should be a function of redshift z (e.g., Tegmark & Peebles 1998; Blanton ct al. 2000; 
Somervillc ct al. 2001). Then, the method proposed here really measures the weighted 
averaged r over the redshift range up to Zg, the source redshift. 

With r = 1, we expect the right hand side of equation (17) (denoted as r^g) to be close 
to 1 if the binning is fine enough. On the other hand, it is observationally impractical if 
a very fine binning is required. In the following, we study how deviates from 1 with 
different binning in the case of linear and deterministic bias, i.e., r = 1. This gives us an 
evaluation on the intrinsic uncertainty when we use equation (17) to measure the r-factor. 

We assume all source galaxies are located at redshift Zg. The foreground-galaxy distri- 
bution is written as 

Pfi^f) = ~ Zi< Zf < Z2 

Z2 — Zi 

— otherwise (18) 

Equal binning is applied in the calculations, i.e., {z2 — zi) = z^^/A'tin where Nbin is the 
number of bins used. A nonlinear power spectrum in the form of Peacock and Dodds (1996) 
is adopted. 

Figure 1 presents the integrand of < Ml^{9c) > [denoted as map{z), see eq. (9) but with 
the integrated variable changed from w to z] with -2s = 1 versus redshift z. Figure la is for a 
A-dominated cold dark matter (ACDM) model with — 0.3, = 0.7, Hq — 67 km/s/Mpc, 
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the shape parameter F = 0.2, and as — 0.93. Different curves correspond to 9c — 1', 3', 5', 

and 10', respectively. Figure lb is for the open CDM model (OCDM) with = 0.3, = 0, 
Hq = 67 km/s/Mpc, F = 0.2, and cg = 0.87. For comparison, we plot in Figure Ic the results 
for both models with 6^ = 5'. It can be seen that the contributions to < M^p{6c) > spread 
over a fairly large range of z. For OCDM, the low redshift contributions are more prominent 
than those of the ACDM model. 

In Figure 2 we show 1 — r^^ (assuming linear and deterministic bias) with Zg = 1 versus 
6c for the ACDM model with Ni,in = 5, 10, and 20, which corresponds to Az = z^ — Z\ = 
0.2, 0.1, and 0.05, respectively. We see that with = 0.2, the deviation of r^^ from 1 is as 
high as about 10%, and 1 — r^^ reaches ~ 7% at 20'. Thus with Nun = 5, the accuracy is not 
good enough to give a sound estimation of r by using equation (17). With N^jn = 10 and 
Az = 0.1, 1 — < 3% for the whole angular range we studied, and at ~ 20', r^^ deviates 
from 1 by only less than 2%. With finer binning of Az — 0.05, 1 — rig < 0.7% for angular 
scales up to 100' with 1 — r^^, < 0.5% around 20'. Therefore, with Az = 0.05, the intrinsic 
uncertainty of equation (17) is tiny, but observationally, it may be difficult to precisely (with 
0.05 accuracy) measure the redshift for a large number of foreground galaxies especially for 
those with redshift close to 1, since most likely only photometric redshifts will be available 
for high-redshift galaxies. Considering both the intrinsic uncertainties and the observational 
realities, it is therefore optimal to use Nbin = 10 {Az = 0.1) in the case of Zg = 1. 

Similar curves for the OCDM model are shown in Figure 3. It is seen that 1 — r^^ is 
larger for the OCDM model than the corresponding 1 — r^^ for the ACDM model. This is 
because as seen in Figure Ic, low-redshift contributions to < M^^ > are more significant 
in the OCDM model than in the ACDM model, and therefore for OCDM, relatively finer 
binning is needed in low redshifts in order to reach the same 1 — r^g as that of ACDM. Still, 
however, with Nj^^ — 10 {Az = 0.1), 1 — r^^ is small enough with 1 — r^^ < 3.5% at 20', and 
1 — rig < 4% for the whole angular range up to 100'. 

In Figure 4 we plot 1 - r^g for the ACDM model {Solid lines) and OCDM model 
{dashed lines) with Zg — 1.5. The two sets of curves correspond, respectively, to Nf,in — 10 
{Az = 0.15) and A^^^ = 15 {Az = 0.1). For N^n = 10 {Az = 0.15), 1 - r^g < 3.5% and 5% 
for the ACDM model and OCDM model, respectively. At 20 arcmin, 1 — rf,^ < 2.5% and 4% 
for the respective two models. For Nbin = 15 {Az = 0.1), the deviations of r^g from 1 are, 
respectively, less than 1.5% and 2% for the two models for the whole angular range consid- 
ered. Thus for Zg = 1.5, Nbin = 15 {Az = 0.1) guarantees that the intrinsic uncertainties in 
terms of r estimation with equation (17) are small. Even with Nbin = 10 (A;^ = 0.15), the 
intrinsic uncertainties are insignificant. 

Figure 5 shows the curves for Zs — 0.5. The results for each of the two models with 
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Nbin — 5 {Az — 0.1) and Nbin — 10 {Az — 0.05) are shown in the plot. For Nbm — 5 and 
Az = 0.1, the deviations are relatively large, with 1— r^^ as high as 8% and 10% for the ACDM 
model and OCDM model, respectively. With Nhin = 10 and Az = 0.05, the deviations are 
less than 2% and 3% for ACDM and OCDM, respectively, with l-r^^ < 1.5% and 2% at 20'. 
Therefore, for Zg = 0.5, it is necessary to go to Az = 0.05 in order to control the intrinsic 
uncertainties in the r-estimation using equation (17). Fortunately, for z < 0.5, the redshifts 
of foreground galaxies can be obtained spectroscopically with very high precision, and it is 
possible to bin them down to Az = 0.05. 

From above analyses, we see that for source galaxies at Zg, in general 10 bins for the 
foreground galaxies would be enough to control the intrinsic uncertainties of equation (17) 
to a reasonably low level. For Zg > 1, the intrinsic uncertainties can be well controlled if the 
bin interval Az — 0.1 is used. 



3. Signal-to-Noise Ratio 

In this section, we discuss the signal-to-noise ratio of r estimated with equation (17) 
considering the intrinsic ellipticity of source galaxies, Poisson noise of the foreground-galaxy 
distribution, and the cosmic variance. 

If we denote 

{[< MaMWec.zf) >]V < N',{e,,zf) >}pf{zf)Azf 

as M^f, then equation (17) becomes 

To estimate the signal-to-noise ratio of r, we assume that the M^^ of different redshift 
bins are independent. This is a reasonable assumption with the bin intervals we considered. 
Further we assume that the dispersion of < Ng{9(., Zf) > can be neglected (van Waerbeke 
1998). For one field, we then have 

^(^'/) ^ 2a{<Ma,{e,)N,{e,,zj) >) 
Ml <Map{9c)Ng{9,,Zf)> ' ^ ^ 

and 

^'(S.,My = S,,a2(My . (21) 
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If E2^(M|^) and < M^p(6c) > are independent, then 

However, it is likely that I^Zf{M^f) and < M^p{9c) > are not independent, so then 
(Taylor 1982) 

In the following, both equation (22) and (Tjnax defined in equation (23) are calculated. 
We use the estimation of van Waerbeke (1998) to calculate a[< Map{9c)Ng(9c, Zf) >], 



I.e., 



a{< MaM)NM. ^f) >) = [ < Ml{e,) > [ < N'^{0,, zf) > "\ (24) 

where Nf, and Nf are the numbers of source and foreground galaxies in one field, respectively, 
(Te represents the intrinsic ellipticity of source galaxies, G = ttO"^ J d'^6Q'^{9,9c) = 6/5 if the 
compensated filter defined previously with / = 1 is used, and G = Tidl J d^9U^{9,9c) = G. 
Within each pair of square brackets of the right-hand side of equation (24), the first term 
represents the uncertainty due to the cosmic variance. For the aperture mass, the intrinsic 
ellipticity of background galaxies contributes to the uncertainty through the second term 
of the first square bracket. The second term in the second square bracket represents the 
statistical uncertainty due to the finite number of foreground galaxies. 

For ct(< >), we have (Schneider et al. 1998) 

_^,<Ml>+^-^ + V2<Ml>)\ , (25) 



r ,.0.. . G 



1/2 



a{< >) ^ 



where //4 is the kurtosis of Map with jiA —< M^p > /(< M^p >)2 — 3, which is if Map is 
Gaussian. On small scales, however, non-linear effects are strong, and Map is non-Gaussian. 

As shown in van Waerbeke (2002), < M^^ > / < > Gaussian can reach about 7 for 9c < 3' 
for ACDM model. On the other hand, on small scales, we expect statistical uncertainties to 
dominate over the cosmic variance, thus, neglecting non-Gaussianity should not affect the 
error estimation significantly. In the following, we ignore the /i4 term in equation (25). 

Figure 6 shows a{r-'^)/r'^ of one field with Zs = 1, A^Mn = 10 {Az = 0.1), Uf, — 
30 gal/ arcmin^, — 0.2 and nf — 5 gal/ arcmin^, where ni, and Uf are the surface number 
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densities of background and foreground galaxies, respectively. The solid lines are for the 
ACDM model, and the dashed lines are for the OCDM model. In each model, the upper 
line corresponds to the result of (y„iax in equation (23) and the lower line corresponds to 
that of equation (22). It is seen that at 9c < 3', the statistical uncertainties arc dominant, 
and oir^^jr^ is larger for the ACDM model. At angular scales Qc > 3', the two models have 
about the same cr{r'^)/r'^, and it approaches a constant when 6c > 10'. 

In Figure 7 we plot a{r'^)/r'^ for a survey of 25 deg^. The meanings of different lines are 
the same as those in Figure 6. We see that o'(r^)/r^ < 20% or a{r)/r = l/2a{r'^)/r^ < 10% 
for 6c < 10'. In other words, the deviation of r from 1 is detectable if it is larger than 10% 
at 6c < 10'. Going deep in a survey, the number density of source galaxies increases and the 
statistical uncertainties decrease. As we have seen, however, on relatively large scales, the 
cosmic variance dominates, and thus a deeper survey will not help to increase the precision 
on the determination of r on those scales. On the other hand, by increasing the survey area, 
(T(r^)/r^ decreases by a same factor over the whole angular range. In Figure 8, we compare 
the effects on (T(r^)/r^ with a deeper survey and with a wider survey for the ACDM model. 
The solid line shows the result of equation (22) with Ub = 30 gal/ arcmin^ and a survey area 
of 25 deg^. The dashed line is for = 60 gal/ arcmin^ and a survey area of 25 deg^. The 
dot-dashed line is for = 30 gal/ arcmin^ and a survey area of 100 deg^. It is clearly seen 
that in terms of determining r, it is more effective to go wide than to go deep in a lensing 
survey. 

Figure 9 shows the results for Zg = 1.5 and Nfym = 15 {Az = 0.1). The other parameters 
are the same as in Figure 7. At small angular scales, {6c < 4'), a{r^)/r^ is smaller in the 
case of = 1.5 because of relatively large cosmic-shear signals compared Zg — 1. A 10% 
detection is reachable for 6c<10' with a survey of 25 deg^. 

In the case of Zg = 0.5, since the shear signals arc very weak, the statistical uncertainties 
are relatively very large. Figure 10 plots cr(r^)/r^ for rib = 30 gal/ arcmin^, Nbm = 10 
{Az — 0.05) and a survey area of 25 deg^ for the two cosmological models. The number 
density of foreground galaxies is Uf — 5 gal/ arcmin^. We see that over the whole angular 
range, a{r'^)/r'^ > 20%. The good strategy is then to do shallow, wide surveys. For the Sloan 
Digital Sky Survey (SDSS), the survey area can be as large as 10^ deg^. The surface number 
density of source galaxies is about a few (e.g., Mckay et al. 2001). In Figure 11 we show 
(T(r^)/r^ for n;, = 5 gal/ arcmin^ and a survey area of 10^ deg^. All the other parameters 
are the same as in Figure 10. It is seen that with such a large survey area, a 10% detection 
of r is possible for all the angular scales we considered (1' — 100'). 

We would like to emphasize that in this section we only aimed at giving an estimation 
on the signal-to-noise ratio of r, i.e., to what level the deviation of r from 1 is detectable. For 
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that, we only calculated the standard deviation of r. Since the statistics of r are certainly 
non-Gaussian, and rather complicated, strictly speaking we cannot draw conclusions on the 
significance level of a measured r based on Gaussian statistics. For example, if we get a 2a 
measurement on r, the Gaussian probability that r is consistent with r = 1 is 5%. For real 
r, the probability can be very different depending on the statistics of r. Nonetheless, we 
expect that our signal-to-noise ratio should give a reasonable estimation of the detectability 
of r. Detailed statistical properties of r can be best studied through numerical simulations, 
which is our planned next step. In addition, in our estimation, we implicitly assumed that 
the uncertainties were small, so that the error propagation rules were applicable. From our 
calculations discussed above, the uncertainties are indeed small on the angular range of less 
than < 10' if the survey area is large enough. 



4. Discussion 

In this paper, we studied the feasibility of extracting information on r from the cross- 
correlation between the foreground-galaxy distribution and cosmic gravitational lensing ef- 
fects under the availability of the redshift information for each foreground galaxy. We showed 
that for ^5 ~ 1, the intrinsic uncertainty in the r-estimation can be well controlled if the fore- 
ground galaxies are binned with A2; = 0.1. This binning requirement is within the accuracy 
range of photometrically determined redshifts, and is thus achievable. For shallow surveys 
with Zs ~ 0.5, a fine binning with Az = 0.05 is necessary. This fine binning is possible, since 
the redshifts of low-redshift galaxies can be obtained spectroscopically. 

We further found that for a moderately deep lensing survey with ~ 1 and = 
30 gal/ arcmin^, a 10% detection on |r — 1| is possible in the angular range of 1' — 10' 
with a survey area of 25 deg^. A deeper survey with Ub = 60 gal/ arcmin^ helps to reduce 
uncertainties only on small angular scales where the noise due to the intrinsic ellipticity of 
source galaxies dominates. The effective way to reduce the uncertainties over a large angular 
range is to increase the survey area. For example, with a survey area of 100 deg^, a 10% 
detection of |r — 1| can be made up to ~ 20'. 

For Zs — 0.5, because of the relatively low cosmic shear signal, the statistical noise from 
the intrinsic ellipticity of source galaxies is much more dominant than in the case of 2;^ ~ 1. 
Thus a very large survey area is needed in order to detect the deviation of r from 1. SDSS 
will eventually have a survey area of about 10^ deg^, and will be very suitable to be used to 
measure r with the method we propose here. 

For the ACDM model, the lensing effects on source galaxies with Zs — 1 mainly come 
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from the matter distribution around redshift z ~ 0.18 — 0.68, with a peak located at z ~ 0.4 

for Qf. = v. The three-dimensional wave- vector k corresponding to ^ ~ 0.4 and 6c = 
V is k ^ 12( Mpch~^)~^. For 6c = 10', the contributions to the lensing effects peak at 
z ~ 0.29, with significant amounts from z ~ 0.1 — 0.57. The corresponding peak scale is 
k ~ 1.7( Mpch^^)^^ . Thus, measuring r in the angular scale range 1' — 10' with Zg = 1, wc 
probe the bias of structure distribution on scales of about 0.05 Mpc/i^^ to about 1 Mpc/i~^ at 
redshift z ~ 0.1 — 0.7. For Zg = 0.5, the lensing effects are mainly from matter fluctuations in 
the redshift range z ~ 0.08 — 0.4 for 6c — l' — 10'. The corresponding scale is ~ 0.06 Mpch"^ 
to ~ 0.6 Mpch-^. 

In the above studies, we have assumed that source galaxies arc all located at the same 
redshift Zg. For future surveys, it is possible to bin source galaxies according to their redshifts, 
and thus our analyses here will be directly applicable. Currently, most surveys select source 
galaxies according to their luminosities. Their redshift distribution is approximately in the 
form Pb{z) oc z'^exp[—{z/zQ)'^] with (3 ~ 1.5. In this case, in order to control well the intrinsic 
uncertainties in the r-cstimation, the redshift range of foreground galaxies has to be extended 
to about < z >= l.bzQ. With zo — 0.7, < z >^ 1. 

With the fast growth of gravitational lensing surveys, we expect, with the additional 
information on the galaxy distribution, that we can obtain a great deal of knowledge on the 
process of galaxy formation, which will further help us to extract cosmological information 
from other types of surveys. 

We thank the referee for the constructive comments that helped us further clarify our 
studies. This research was supported by the National Science Foundation of China, under 
Grant 10243006. 
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Fig. 1. — Plot of map{z)^ the integrand of < M^^ >, vs. redshift z. (a) For tlie ACDM 
model. The sohd line is for — 1', the dotted line is for — 3', the dashed hne is for 
= 5', and the dot-dashed hne is for = 10'. (b) For the OCDM model, (c) For = 5', 
where the solid hne is for the ACDM model, and the dashed hne is for the OCDM model. 
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Fig. 2. — Intrinsic uncertainty 1 — r^^ vs. the angular scale 6c (arcmin) with Zg = 1 for the 
ACDM model. The solid line is for Az = 0.05, the dotted line is for Az = 0.1, and the 
dashed line is for Az = 0.2. 
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Fig. 3. — Intrinsic uncertainty 1 — vs. the angular scale 9^, (arcmin) with -z^ = 1 for the 
OCDM model. The solid line is for = 0.05, the dotted line is for Az = 0.1, and the 
dashed line is for Az = 0.2. 
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Fig. 4. — Intrinsic uncertainty 1 — r\g vs. the angular scale Qc (arcmin) with Zg — 1.5. The 
solid lines are for the ACDM model, and the dashed lines are for the OCDM model. In each 
pair of lines, the upper one is for iS.z — 0.15, and the lower one is for iS.z — 0.1. 



= ACDK/ 


1 1 1 1 1 

1, OCDM 


1 1 1 1 
, z =1 5 


- 




Az = 0.1 5 






" — ■ 


" ^ ' 






; / 
/ 








/ 

r / 

: / 

: / / 
- 1 / 

y ,'-~ 
-_i / y 








Az-0.1 








Az = 

















III, 


1 , , , 1 , 





-20- 



g.5 



00 

CD 



0.12 
0.10 

0.08 

0.06 

0.04 

0.02 

0.00 



ACDI 



- OCDM, Z3=0.5 
Az = 0.1 



Az = 0^ 



- / 



Az = 0.05 
"Az~=0.05 



20 



40 60 



80 



00 



Fig. 5. — Intrinsic uncertainty 1 — r^^ vs. the angular scale 9c (arcmin) with Zs = 0.5. The 
solid lines are for the ACDM model, and the dashed lines are for the OCDM model. In each 
pair of lines, the upper one is for Az — 0.1, and the lower one is for Az — 0.05. 
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Fig. 6. — Plot of a{r'^)/r'^ vs. 6c (arcmin) for one field with Zg = 1, n(, = 30 gal/arcmin^, 
rif = 5 gal/arcmin^, and — 0.2. The sohd hnes are for the ACDM model, and the dashed 
lines are for the OCDM model. In each pair of lines, the upper one corresponds to a^ax in 
eq. (23), and the lower one corresponds to eq. (22). 
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Fig. 7. — Plot of a{r'^)/r'^ vs. 9c (arcmin) with a survey area of 25 deg^. All the other 
parameters are the same as in Fig.6. 
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Fig. 8. — Plot of a{r'^)/r'^ vs. 9c (arcmin) for the ACDM model with 



1 (eq.[22]). 



The sohd hne is for rib — 30 gal/arcmin and the survey area of 25 deg , the dashed hne 
is for Ub — 60 gal/arcmin^ and a survey area of 25 deg^, and the dot-dashed line is for 
rib — 30 gal/arcmin^ and the survey area of 100 deg^. 
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Fig. 9. — Plot of a{r'^)/r'^ vs. 9c (arcmin) with Zg = 1.5, = 30 gal/arcmin , Uf = 
5 gal/arcmin^, and = 0.2. and a survey area of 25 deg^. The meanings of the lines are 
the same as in Fig. 6. 




Fig. 10. — Plot of a{r'^)/r'^ vs. 6c (arcmin) with Zs = 0.5, rib = 30 gal/arcmin^, Uf = 
5 gal/arcmin^, = 0.2, and a survey area of 25 deg^. The meanings of the hnes are the 
same as in Fig. 6. 
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Fig. 11. — Plot of air^^jr^ vs. (arcmin) with Zs = 0.5, nj, = 5 gal/arcmin^, Uf = 
5 gal/arcmin^, = 0.2, and a survey area of 10"^ deg^. The meanings of the lines are the 
same as in Fig. 6. 



